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Abstract 

A family of diffusion-annihilation processes is introduced, which is exactly 
solvable. This family contains parameters that control the diffusion- and 
annihilation- rates. The solution is based on the Bethe ansatz and using 
special boundary conditions to represent the reaction. The processes are 
investigated, both on the lattice and on the continuum. Special cases of 
this family of processes are the simple exclusion process and the drop-push 
model. 
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1 Introduction 



In recent years, the asymmetric exclusion process and the problems related to 
it, including for example the bipolimcrization 0], dynamical models of interface 



shocks JS), |6[, have been extensively studied. The dynamical properties of this 
model have been studied in [6-8]. As the results obtained by approaches like 
mean field are not reliable in one dimension, it is useful to introduce solvable 
models and analytic methods to extract exact physical results. Among these 
methodes is the coordinate Bethe-ansatz, which was used in || to solve the 
asymmetric simple exclusion process on a one-dimensional lattice. In |To[ . a 
similar technique was used to solve the drop-push model fll|| , and a generalized 
one-parameter model interpolating between the asymmetric simple exclusion 
model and the drop-push model. In this family was further generalized 
to a family of processes with arbitrary left- and right- diffusion rates. All of 
these models were lattice models. Finally, the behaviour of latter model on 
continuum was investigated in [T^ | . The continuum models of this kind are also 
investigated in . 

In the generalized model interpolating between the asymmetric simple ex- 
clusion model and the drop-push model |l2|, [k| , there are two paprameters 
A and /x, which control the pushing rate. Normalizing the diffusion rate to 
one, it is seen that the sum of these two parameters should be one to ensure 
the conseravtion of probability. These two parameters appear in the boundary 
condition used instead of the reaction. A question arises that what happens if 
one violates this conservation of probability. This is what is investigated in the 
present paper. 

The scheme of the paper is the following. In section 2, the allowed boundary 
conditions are investigated. It is shown that if A + /i < 1 , then the number of 
the particles will be decreasing, that is, there is an annihilation process as well. 
It is shown that one can in fact write a two - particle to one particle annihilation 
process which results in such a boundary condition. 

In section 3, the Bethe-ansatz solution for the ./V-particle probability of this 
process is obtained, and its large-time behaviour is investigated. This is done 
for the process on the lattice as well as on the continuum. Finally, in section 4 
the special case of the two-particle initial condition is fully investigated, and it 
is explicitely shown that at large times, there remains only one particle. 

2 Boundary conditions 

Consider the following master equation for an asymmetric exclusion process. 




—P(x 1: x 2 , • • • , 



x N ;t) 



P(x 1 - l,x 2 , ■ ■ ■ ,x N ;t) 
+P(xi,x 2 - 1, ■ • ■ ,x N ;t) H 
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+P(xi,x 2 , ■ • • - 1;£) 

-NP(xi,X2,- ■ - ,XN)t). 



(1) 



This equation describes a collection of N particles drifting to the right with unit 
rate. If the particles are to exclude each other, that is if no two particles are to 
occupy the same site, then ([!]) is valid only for 

Xi < x l+ i - 1. (2) 

One can, however, assume that ([j]) is correct for all of the physical region Xi < 
Xi + i, and impose certain boundary conditions for xi = Xi+%. Note that if 
Xi = Xi+i — 1 for some i, then in the right-hand side of (|l|) there will be terms 
with Xi — Xi+i, which is out of the physical region. The boundary condition 
determines the nature of the interaction between particles. But what are the 
allowed boundary conditions? Let's rewrite (|l|) for the case of two particles and 
use the conservation of probability. We arrive at 

]T p (zi>z2;t) 

X 2 2)1<I2 
X 

+ 12 12 P{xi,x 2 ;t) 

X 2 Xx<X 2 
X 

x 2 X\<Xi 

-^P(x,x + l;t)+^P(x,x;t). (3) 

X X 

If the right-hand side of (||) is to be identical to zero, then P(x, x) should be a 
linear combination of P(x — 1 + i, x + i)'s. This may work for the two-particle 
process, but in many-partcles processes it may introduce terms with Xi > a^+i, 
which need additional boundary conditions. The only exception is when P(x, x) 
is a linear combination of only P(x, x + 1) and P(x — l,x). So one can write 

P(x,x)=XP(x,x+l)+fj,P(x-l,x). (4) 

Inserting this in ([|), one arrives at 

§- t 12 12 P(.x 1 ,x 2 ;t) = (\ + ii-l)'22P(x,x+l;t). (5) 

X 2 Xi<X 2 X 

In order that the right-hand side of (||) be zero, one must impose 

A + /! = !. (6) 



dl!2 12 P^u^t) = 

x 2 xi<x 2 
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This is the boundary condition used in 10 [l2, 13 

This kind of boundary condition ensures the conservation of particle number. 
But in a process where annihilation exists as well, the number of the particles 
is not conserved; it is decreasing. It is seen that if 

X + fKl, (7) 

then the probability of finding two particles is decreasing. Suppose one begins 
with two particles. They drift to right with unit rate. If they meet each other, 
either the left particle is stopped, or one of them is annihilated. That is, we 
have the following processes: 

A% -> %A, with rate 1, 
A A -> $A, with rate a, 

AA A®, with rate j3. (8) 

In this case, still the more-than-two particle densitiies are zero, since no par- 
ticles are generated during the process. But the two-particle density does not 
determine the one-particle density. And the summation of the former need not 
be a constant (one). It is seen that the master equation for the two-particle 
probability is 

d 

—P(x 1 ,x 2 ;t) = P(xi - l,x 2 ;t) + P(xi,x 2 - l;t) 

-2P{x 1 ,x 2 ;t), x 1 <x 2 -l, (9) 

and 

d 

—P(x, x + 1; t) = P(x - 1, x + 1; t) - (1 + a + /3)P(x, x + 1; t). (10) 
ot 

But ([n]) is the same as (||), provided one uses the boundary condition 

P(x,x) = XP(x,x + l), (11) 

with 

X=l-(a + P). (12) 

So the difference 1 — A is in fact related to the annihilation rate, as expected. 
There is one other thing to be noted. As the number of particles is not conserved, 
one cannot calculate the one-particle probability by a simple summation of the 
two-particle probability. That is, 

PW^p(i,i,) + ^%4 (13) 

y>x y<.% 

In fact, for the process described, the particles interact and annihilate each 
other, until there remains only one particle. This means that at t — > oo, there 



3 



will be only one particle. So the more-than-one-particle probabilities will tend 
to zero, whereas the summation of the one-particle probability tends to one. 
However, if the initial condition is that there are N particles, one can write dif- 
ferential equations for n-particle probabilities in which n-particle probabilities 
and n + 1-particle probabilities occur (if n < N). For n > N, the n-particle 
probability is identically zero, and the equation for the iV-particle probability 
is closed. So, in principle, one can find the ./V-particle probability first and 
use this to find less-than-iV-particle probabilities. To be specific, the evolution 
equation for the one-particle probability is 

^P(x;t) = P(x-l;t)[l-P(x;t)]-P(x;t)[l-P(x + l;t)] 

-aP(x,x + l;t) - j3P(x - l,x;t). (14) 

3 Bethe— ansatz solution for the N— particle prob- 
ability 

Consider the Master equation (Q), with the boundary condition 

P{---,x,x,---) = XP(---,x,x +!,••■), (15) 



where A < 1. Following g, [1C| |12|, |13|, one can obtain the conditional probability 
using the Bethe ansatz. Writing 

P(x;i)=e £t vl>(x), (16) 

and 

*(x)^ie w ( p )' s , (17) 

a 

where the summation runs over the elements of the permutation group, one 
arrives at 

E = -N + J2 e ~ iPi > ( 18 ) 

3 

and 

A aai = SMpOX^+OIAt, (19) 

where a is that permutation which only interchanges pi and Pi+\ . One also finds 
that 

1 - Xe lPk 
1 - \e % Pi 

This is the same as what found in |l^, [l2| with /.i = 0, and if one puts A = 1, 
the result of H is obtained. The conditional probability is thus written as 



Sjk ■= S(pj,p k ) = 7—-. (20) 



r rfN 

P(x; t|y; 0) = J ^(x)^-^, (21) 
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where \l/ is defined as ( |l7| ) with ^identity = 1- This looks like similar to what 
obtained in || [l0| [l^]. There is, however, a difference. As A < 1, there is no 
pole in S, and hence in A. So for large times, when the probability distribution 
becomes smooth and its Fourier-transform for large frequencies tends to vanish, 
one can put Pj = in S as an approximation to arrive at 

S«-l, (22) 

and 

A a » {~iy. (23) 
One can also approximate -E(x) as 



^(p)«E^-yJ- ( 24 ) 

So, for large times, 

P(x;i|y;0) « ^-L_ ^(-l)^-E,[%--fe)-*] 2 /(2*)_ (2g) 

It is clearly seen that the integral of this distribution over the physical region 
tends to zero as t — > oo. This should be the case, since the number of the 
particles does not remain constant and decreases. 

Using the boundary condition (^) with ((?]), doesn't change the results dras- 
tically. In fact, E doesn't change at all, while S is changed to 

Sjk = - ^— . (26) 

The approximate result for large times does not depend on A or /i, so long as 
their sum is less than 1. 

One can also investigate the continuous-space form of the evolution. Fol- 
lowing p3[ , the master equation is changed to 

3 3 

and the boundary condition to 

(1 - A - n - Xd j+1 + iidj)P\ Xj+1=Xj = 0. (28) 

Using the Galilean transformation xi — > Xi + vt and t—*t, the master equation 
( p7| ) is simplified to 

|p(x;i) = iv 2 P(x;t), (29) 



2 
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E=- 

2 

j 



Using a Bethe-ansatz solution like (|l^) and ([j~7|), one arrives at 

3 

1 - A - \i - ^Ap fc + z/xpj 
1 — A — /i — iApj + i/ij»fe 

For large times, one can approximate S to — 1, and arrive at a result similar to 
(p5|). The difference is that in the exponent the term xj — t — a(yj) is replaced 
by Xj — cr(yj), as the Galilean transformation used has canceled the drift from 
the master equation. 



and 

Sjk 



4 Two— particle system and the exact solution 



As it was seen in the previous section, the conditional probability for the two- 
particle system described by (||) and (|TT|) is 



P(x;t|y;0) = 1 



1 - Xe^ 



,i(pix 2 +P2Xi) 



(32) 



1 - Ae*Pi 

where E is obtained from (|18|). This integration is easily done and the result is 



P(x;t|y;0) = e 



Oi - yiY- {x2 - 2/2)! 

-2/ 



2 ^ ^— -A' [ —1 



1=0 



(l + x 2 - J/i)! (xi - y 2 )! 



Af 



xi - 2/2 + 1 



133) 



Another interesting quantity is the average number of the particles. This is 
equal to the summation of the one-particle probability: 



N(t) :=53P(x;i). 

X 

Using (pT|), one arrives at 

N = -(a + f3)J2 P(x - l,x;t) 



(34) 



1 - A 



(35) 



G 



The right-hand side can be calculated using the Bethe-ansatz solution directly. 
Using ([32|), one has 



J 2ir 1 - Ae*P 

2n 



e %{m-l)p _ gj(m+l)p 



xE^ r 

m=0 
00 

A e 2t A m [I J , 2 _ 3/1 _|_ m _i(2i) — L y2 _j, 1+m+ i(2t)] 



A £ e"-A 



m=0 



^y2-yi+m('2t) 1 



(36) 



where I denotes the modified Besscl function. This can be inserted in ( pq ) to 
obtain 

JV(t) - jv( )-i^ f dt'Y^p&x-j) 

00 _ t 

= JV(0) - (1 - A) £ / A 
This is simplified for t — ► 00. Using 



3/2— Ki+m 



(2f).(37) 



e~ s 1 

ds I n (s) = -, 

s 71 



(38) 



one arrives at 



7V(co) = N(0) - (1 - A) 
= N(0)-1. 



m— 



But note that 



7V(0) = 2, 

since at the beginning there were two particles at y\ and j/2- That is, 
P(x;0)=8{x-yi)+6(x-y 2 ). 



(39) 
(40) 

(41) 



So, at t — > 00, there remains only one particle, as one of the two particles has 
been annihilated. 



The continuous-space analogue of this model can also be solved easily. Using 
9|) as the master equation, and ( p0| ) and ( |3l| ) (with /i = 0), one is led to 



P(x;t|y;0) 



_4tt 2 

e 



i(pix 1 +p 2 x 2 ) _ 1 A i^P2 ( J,(p 1 x 2 +P2Xl) 

1 — A — iXpi 



(42) 



Using the change of variable p := p\ + i(l — A)/A in the second integral, ( p2| ) is 
written as 

P(x:t|y;0) = -^ e -[(^-w) 2 +(^-^) 2 l/( 2t ) 



7^ ^ 

/8?rf 

X 1 -1 + erf 



e -(Kl-y2)V(2t) e [2^(2:2-ai)+tA- ! ]/2 



—?={X2 - yi+tA) 
V 2l 



where 



.4 



1 - A 
A ' 



(43) 



(44) 



One notes that at t — ► 00, the conditional probability is simplified. We have 



1 — erf(x) 



for large a; 



(45) 



From this, it is seen that at t — ► 00 



P(x;t|y;0)«^{ e -^-^ 



i-+(^2-y2) 2 ]/(2t) _ e -[( :ci - y2 ) 2 + ( a;2 - yi ) 2 ]/(2t)| _ 



(46) 

This is a special case of what obtained in the previous section. 

Another quantity to be considered is the one-point probability. In the con- 
tinuum limit, and after performing the Galilean transformation, ( |l4| ) becomes 



From this, using the boundary condition (|28|), with /i = 0, one arrives at 

1 - A 



tf(t) = - 



A 



dx P(x, x; t), 



where 



N(t) := / dx P{x-t). 



(48) 
(49) 
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Using (^2|), the integral at the right-hand side is calculated to be 



To obtain N(po), one integrates ( |50| ) from to oo. This results in 

f A f dv e ip ( V2 ~ Vl ^ 

dt / dx P(x,x;t) = -P ' ' 



7r J ip 1 — A — iXp 
A 



1-A 



(51) 



The symbol P denotes the Cauchy's principle value, and use has been made of 
the fact that yi > y\. From this, it is found that 

7V(oo) = AT(0) - 1 

= 1- (52) 

This is the same result obtained for the lattice, as expected. 
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